International Journal of Power Electronics and Drive Systems (IJPEDS) 

Vol. 9, No. 3, September 2018, pp. 1321 - 1329 
ISSN: 2088-8694, DOI: 10.11591/ijpds.v9.i3.ppl321-1329 


□ 


1321 


A New Adaptive Anti-Windup Controller for Wind Energy 
Conversion System Based on PMSG 

Ed-dahmani Chaflk, Mahmoudi Hassane, Bakouri Anass, and El Azzaoui Marouane 

Power Electronics and Control Laboratory, Electric Department, Mohammed V University, Morocco 


Article Info 


ABSTRACT 


In this paper, an adaptive anti-windup control strategy for permanent magnet synchronous 
generator dedicated for wind energy conversion systems. The proposed control has the ad¬ 
vantage to suppress the performance deterioration caused by the overshooting phenomenon, 
and optimize the controller gains using the particle swarm optimization algorithm. The 
scheme of the speed controller is implemented on field orientation control in the generator 
side converter. A simulation of the proposed scheme is carried out in SIMULINK-MATLAB 
in order to evaluate the effectiveness of the control against the saturation and the parameter 
optimization. 
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1. INTRODUCTION 

In the last decades, there has been a growing interest in wind turbines. Electrical generators and control 
strategies should respond to the needs of wind power applications. The permanent magnet synchronous generator 
(PMSG) and doubly fed induction generator (DEIG) are widely used in wind energy conversion systems (WECS) 
because they offer the possibility to work with variable wind speed [1]. Also, they can offer an improvement on 
production of wind energy, and the ability to achieve maximum energy conversion efficiency. 

Comparing to doubly fed induction generator (DEIG), the PMSG can provide a high-efficiency and high 
reliability power generation, low maintenance required, and the electrical losses in the rotor are eliminated [2],[3]. 
Due to the mentioned advantages, the PMSG becomes an interesting solution for wind turbine applications. Eor 
a special architecture with a high number of poles pairs, the PMSG offers the possibility to eliminate the gearbox 
system as shown in Eigure 1, that allows a higher efficiency. Several authors [4],[5] have discussed about the Variable 
wind speed conversion systems based on PMSG. 

The conventional control strategy of PMSG is based on field oriented control (EOC) with PI controller [6]. 
This controller is easy to apply, and present a good performance in linear region [7], [8]. But, it suffers from non-linear 
effects such as saturation, when the saturation is neglected in the design phase it causes instability during closed-loop 
where the output system response takes a long time to stabilize in the steady state. In other words, this phenomenon is 
caused by the windup integrator contained in the PI controller, which keeps integrating the tracking error even if the 
input is saturating [9]-[10]. 

In order to overcome the windup phenomenon, several researches have proposed anti-windup techniques to 
deal with input saturation, where some proposals involve a complicated design [11],[12]. Thus, the average strategy 
to handle the integrator windup by tuning the controller disregarding the saturation caused by the integrator, and then 
adds an anti-windup compensator to avoid performance graduation [13]. Basically, the classical anti-windup strategies 
come with two different approaches, namely, conditional integration and tracking back-calculation. 

In this paper, the speed control is based on adaptive anti-windup PI controller. The adaptation controller 
gains are gotten with search technique known as Particle Swarm Optimization technique (PSO), this technique was 
developed by Eberhart and Kennedy in 1995 [14],[15]. This technique is an optimization tool based on population. 
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and the system is initialized with a population of random solutions and can search for optima by the updating of 
generations. The PSO algorithm has been used in power system for tuning control purpose in [16],[17],[18]. 

The aim of this paper is to implement, discuss and compare the racking performance of adaptive anti-windup 
speed controller with conventional linear PI. The proposed scheme of anti-windup eliminates the overshooting with 
a simple structure existing in PI controllers, and guarantees the independence between the desired performance and 
operating conditions. The simulation is realized with Simulink in order to verify the performance improvement com¬ 
paring to conventional PI controller. 


2. MODELING THE WECS CONCEPTS 
2.1. Wind turbine model 

A wind turbine is composed of many parts to convert a kinetic-to-electrical energy. The mechanical power 
and torque delivered by a wind turbine is given by (1) and (2). 

P^ = O.bpCpTTR^V^ (1) 

Tt = 0.5pCpTrR^— (2) 

Wt 

Cp = 0.5 — 0.4/3 — 5^ • exp ^ ^ ^ + 0.0068A (3) 


With- — — _ - _ 0035 

Ai “ A+O.08/3 

p is the air density, R is the blade radius, v and Wt are respectively the wind and turbine speed, Cp is the power 
coefficient, Pt is the turbine power, and Tt is the turbine torque. The tip speed ratio (TSR: A) is an important parameter 
in wind energy systems. It is defined as the ratio of the blade speed to the speed of incoming wind. 




Figure 1. WECS based on direct drive PMSG. 


Figure 2. Turbine power evolution curve. 


2.2. PMSG Model 


Similarly to the induction generator (IG), the construction of the stator in PMSG is essentially the same. On 
the other side, the rotor magnetic flux is constant and generated by permanent magnets [19]. Also, depending on 
magnets architecture, the PMSG can be classified into surface-mounted and inset PM generators. The mathematical 
model of PMSG in synchronous rotating dq-reference is given by the following equations. 


Vsd = 


^s-Cd Cd “f WrLq.isq 

Rs-isq ^rCd-’i’sd T 

3P 

7e ~ Cq^igd-Cq “t” •’^sq) 


jdWrr. 

dt 


= Tt-Te- BWr, 


(4) 

(5) 

(6) 


Where, Rg is the stator resistance, Ldq are the inductances in dq-reference, is the PM-flux, P is the pole 
pair number, and Wm are electrical and mechanical speed, VdqJdq are the stator voltages and currents components 
in dq-reference, Te and Tt are respectively the electromagnetic and turbine torque, J and B are respectively the 
equivalent system inertia and viscous damping. 
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The equations (4) represent the electrical behavior of the PMSG which gives the possibility to control the 
dq-currents components. Also, the speed controller is based on (6). It should be noted that the electromagnetic torque 
Te may be controlled directly by the quadrature current component in case of surface mounted PM or zero direct 
current control (ZDC). Then, the expression of Tg is given by (7). 

3P 

Ze — 4^ f (7) 


3. WECS CONTROL STRATEGIES 
3.1. Wind turbine control 

The purpose by turbine control is to produce a maximum power, this is achieved for a particular value of 
power coefficient, called Cp-max, the maximal power coefficient is set to 0.48, and it is given for a specified TSR 
Aopi =8.1 with pitch angle /3 = 0. The maximum power operation can be achieved with optimal torque control 
according to (9). 


^opt - 


R-'^m—opt 


Rt — opt — ^-^pCp—max 




^opt 


( 8 ) 

(9) 


Also, to guarantee the safety working of the turbine against strong wind, a pitch angle controller is imple¬ 
mented in Figure 4. The working principle of this controller is by comparing the generated power Pt with the nominal 
power, then it generates a pitch angle response. 



Shaft 

Figure 3. Turbine model diagram scheme. 
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I W 

controller 



Figure 4. pitch angle controller scheme 


3.2. Speed control of PMSG 

The windup phenomenon comes when the output current command of the speed controller is limited to a 
maximum value, and the integral component becomes very large because is not compatible with the plant input. 
Which causes a large overshoot and slow settling time in the speed response. 

In this paper, the speed-loop controller has been designed using an adaptive gains anti-windup algorithm. It has 
an advantage to eliminate the overshooting phenomenon during operation that can deteriorate the power generation 
performance. The aim of anti-windup control for nonlinear system with saturating actuators is to modify the control 
in order to limit the overshoot. 
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The proposed anti-windup strategy by [7] as shown in fig.5 can switch between the P and PI modes according 
to the working states. For the PI mode, it is necessary to initialize the integrator, and in case of P mode the initial value 
of is5_i(0) is inserted, where the PI mode can utilize this value. A low pass filter is used to avoid an abrupt change of 
current by loading the initial value isq-i{0) in the LPF. 


3.3. Closed loop identification 

According to the mechanical equations of the PMSG (6),(7), and for simplicity, the viscous damping B is 
neglected, the new expression of mechanical equation of WECS is given as follow 


dWm 

=Kt.is„-Tt 

dt t sq t 


( 10 ) 


The quadrature stator current component {igq) can be written for the PI controller as (12). 


{ '^sq — ^sq—i '^sq—p 

isq—p — d^p i^'^m—ref ^va) ( 11 ) 

i — — (w ( — W ) + 

Where, Wm-ref is the reference speed of the PMSG which is the optimal speed derived from the MPPT of 
turbine, igq-p and igq-i are the proportional and integral components of igq, Kp and Ki are the PI controller gains. 
By implementing the Laplace transform on (10) and substituting by (11), the transfer function of the closed loop for 
the mechanical is expressed by (12). 


J {sWm - Wm{0)) 





(w-m—re f '^m) T 


tsg(O) 

S 


( 12 ) 


Wm(0) and isq(O) denotes the initial states values of mechanical speed and quadrature stator current. Then, 
the mechanical speed Wm(s) can expressed as a function of the input arguments shown in (13). 



r Ki.Kt+sKt.Kj, n 

T 



Js^-\-KtKp.s-\-KtKi 


Wm—ref 
'i'sq—i ( 0 ) 
( 0 ) 
Tt 

Wm{s) = 

Js-^+KtK^.s+KtKi 

Js 



Js-^+KtKp.s+KtKi 



Ijs’^+KtKp.s+KtKiJ 




(13) 


In the steady state, according to (10) and (11), the proportional and integral terms of quadrature current are 
expressed as: 

{ ^sq—psteady 0 
'^sq—isteady — 


(14) 


When the PI mode is activated, the initial integral term of igq is defined by (15). 


Tt 

^sq—i(0) = TTi (^W'm—ref (15) 

The first term of (15) signify the required current at steady state according to (14), for the second term it is 
assigned to the compensation term for the overshoot. With m is the anti-windup controller gains. Substituting (15) in 

(13). 


Kt{Ki-\-s{Kp — ra)) Kt.m-\-J.s 


W^rn—ref 

Js^ + KtKp.s+KtKi Js'^+KtKp.s+KtKi _ 


_ ^m(O) _ 


Then, (16) can be simplified during the PI mode to the final expression as: 


Wm{s) - 


Wm(0) 


Kt{Ki + s{Kp- m)) / ^ ■lUm(O) 

Js2 + KtKp.s + KtK, ■ s 


(16) 


(17) 
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3.4. Controller design 

In order to establish a good performance for speed controller, the anti-windup controller gain m should be 
determined in the PI mode. According to (12), the transfer function can be rewritten as 

. Kt{Ki + s{Kp-m)) ^ pi.p2 s-z 

Js"^ + KtKp.s + KfKi z ' {s - pi){s - P2) 

Where z is the zero, and pi 2 are the poles of transfer function. 


z = — 


Ki 

K„-r 


-Kt.Kj,±J{Kt.Kp)^-4J.Kt.Ki 

Pi,2 = —-------- 


2J 


(19) 




Figure 5. Diagram scheme of the proposed anti-windup Figure 6. Global structure of the generator side based on 
speed controller. the proposed control for the mechanical speed. 

As mentioned in [7] and [9], the pole is a function of PI gains, and the zero-location determined by the anti¬ 
windup and PI gains. In order to simplify the transfer function to hrst order as shown in (21), it is assumed that z = pi 
in such away \pi\ < |p 2 1, which represents a hrst order low pass hlter without saturation. In that case, the anti-windup 
gain is given by (20). If the gain is smaller than the specihed value, a higher overshoot is provided. Flse, a large value 
can result a slow response. 


m = Kp -\ - - 

Pi 

(20) 

H{s) = 

(21) 


P2- S 


The next step is to determine the initial value of the input arguments. In transition state between the P mode 
to PI mode, the initial integral term of q-axis current component is given as follow: 


i(0) — ^max {^^m—ref tU7n(0)) 


( 22 ) 


Where Imax is the maximum limited current. Also, according to (10), the turbine torque can be expressed by 
in steady state. Using (15) and (22), the new expression of jsq_i(0) is given by 


isg— 2 ( 0 ) — 


^p-^sq—isteady TTl.Ijy 

Kp — m 


(23) 


When the P mode is selected, the initial integral current isq-i{0) is generated by the low-pass hlter. The 
initial mechanical speed Wm(0) can be expressed as follow 


tU7n(0) — 'dim—ref T 


'^sq—isteady “b ^max 


Kp — m 


(24) 


3.5. Particle Swarm Optimization 

In the previous part, the controller has been designed. The anti-windup controller parameters gains, as Kp, 
Ki and m determine the performance system. The selection of parameters is a task that can be difficult. In order to 
guarantee a fast-dynamic response and optimal control, PSO comes to solve the problem of parameters estimation. 
PSO derived from social, psychological theory. It imitates the natural process of group communication to share 
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individual experience flocking, migrating, or hunting. Basically, it searches for the optimal solution from a population 
of moving particles. In PSO, starting with a randomly initialized population called a swarm, each member called 
particle flies through the searching space, where is positioned by Xi vector, evaluating the fitness, and remember the 
best position Xgbest on which it has the best fitness. This information is shared by all particles and adjust their positions 
Xi, and velocities Vi according to the information. The velocity adjustment is based on the historical behaviors of the 
particles themselves and their companions. The particles tend to fly better to best the positions [17],[18]. The velocity 
and current position respectively of every particle are evaluated by (25) and (26). 


Vi{t + 1) = W.Vi{t) + a [ri [Xpbest - Xi{t)) + r 2 {Xgbest - X^{t))] 


(25) 


+ 1) = Xi{t) + Vi{t) 


(26) 


Where, t is time step, w the inertia weight factor, a acceleration constant, ri, r 2 are random functions in the 
range of [0,1], Xi the position of particle, Xptest the best previous position of particle, Xgbest the position of 
best particle among the entire population, and Vi the velocity for the particle. 

The adaptive weighted PSO has been proposed in (27) to improve the reaching capability. 




(27) 


With Nt indicate the iterations number, t is the current step, and oq is a constant in [0.5,1]. It should be noted 
that the inertia weight changes at every step by (28). 


w = Wo + fa (1 - Wo) 


(28) 


With, Wo is a positive constant chosen in [0.5,1], and is a random function in the range of [0,1]. 


4. SIMULATION RESULTS 

To verify the effectiveness of the PI anti-windup speed controller, a simulation of the proposed scheme is car¬ 
ried out in Simulink. The PMSG and turbine data are listed in Table 1 .In order to evaluate the controller performance 
in extreme cases, a step change is applied for the speed reference and load torque. Figure 7 shows a comparison of the 
tracking performance of mechanical speed responses by the conventional PI and anti-windup controllers with different 
iteration numbers for a step changing speed response 70rad/s —)■ 157rad/s —)■ 120rad/s. By comparing the speed 
responses, in the adaptive anti-windup controller, the saturation input is limited which can guarantee a better stability 
and high tracking performance. Also, compared to the conventional PI controller, the PSO improve the controller 
performance, by selecting the optimal gains. With the proposed method, the steady-state is quickly established for 
an optimal value of the iterations number.Therefore, the anti-windup controller provides the optimal dynamic perfor¬ 
mance in term of convergence, saturation, and robustness compared to conventional PI controller. The simulation of 
generator side converter for the WECS is based on the wind speed profile of Figure 8, it should be noted that the 
nominal speed of the wind is chosen Vnom = 12m.s~^. From Figure 9(a)-9(b), the power coefficient and tip speed 
ratio-TSR are maintained at their optimal values by the MPPT control. The pitch angle controller is activated when 
the wind speed exceeds the nominal speed. Then, the power coefficient and tip speed ratio are decreased in order to 
kept the extracted power at the nominal value. Figure 9(c) represents the pitch angle controller response for variable 
wind speed. The mechanical turbine power is illustrated in Figure 9(d) changed according to the wind speed variation, 
also it is maintained in nominal state when the pitch angle controller is activated. 
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Figure 7. Comparison Simulation responses for conven- Figure 8. wind speed profile, 

tional PI speed-controller and the speed anti-windup con¬ 
troller with different iterations number Nt 




Time(s) Time(s) 


(a) Power coefficient response 


(b) Tip speed ratio response 




Time(s) Time{s) 

(c) Pitch angle response (d) Turbine power response 


Figure 9. Turbine dynamic performance using the MPPT and pitch angle controllers 


Table 1. Turbine and PMSG data 

Nominal power 

Turbine radius 

Air density 

Gearbox gain 

Equivalent system inertia 

Maximum power coefficient 

Optimal speed ratio 

Nominal current 

Stator resistance 

Stator inductance 

PM flux 

Pole pairs number 
Nominal speed 
Nominal frequency 


Pn 

l.7kW 

R 

1.04m 

P 

\.22kg/mf 

G 

1.7 

J 

0.35N .m.rad~^ .s'^ 

''p—max 

0.48 

^opt 

8.1 

^nom 

5A 

Rs 

2.7CI 

Ld,q 

3AmH 

(j)f 

OMlWb 

P 

4 

^nom 

157.lrad.s~^ 

fr 

lOORz 


Figure 10(a) shows the high tracking performance of mechanical speed response of the proposed controller 
under a variable turbine speed, it should be noted that the refence speed is given by the MPPT bloc controller. Also, the 
mechanical speed response is stable and tracks the reference value by using the selected value of the PSO algorithm. 
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In Figure 10(b), the electromagnetic torque is identical to the reference value. When, the pitch control is activated 
the PMSG speed and electromagnetic torque are kept at the nominal values, which implies that the extracted power 
is maximal. Figure 10(c) shows current responses in dq frame. The FOC with ZDC is applied, where the reference 
current component of d-axis is set to zero {isd-ref = 0) as shown in Figure 6, and the quadrature current is propor¬ 
tional to the turbine torque as mentioned in (7). For the three-phase stator current response is shown in Figure 10(d). 
Where, the current amplitude and frequency are proportionals respectively to electromagnetic torque Tg and generator 
velocity Wm- 



Time (s) 


(a) Mechanical speed response of PMSG 



10 



0 - 

-2 1 ^^^^^^^^^- 

0 2 4 6 8 10 12 14 16 18 20 

Time (s) 



(c) Current response in dq plan 


(d) Stator current response 


Figure 10. PMSG dynamic performance with the anti-windup speed controller 


5. CONCLUSION 

The adaptive anti-windup was proposed in this paper to replace the conventional PI controller for the speed 
controller in the direct drive PMSG held oriented control. The initial values of the integrator current and mechanical 
speed are determined. The PSO algorithm is used to estimate the optimal parameters of the proposed controller, which 
gives a high tracking and dynamic performance, and fast response with least overshoot. The anti-windup controller 
is designed to give a best tracking speed comparing to the conventional linear PI controllers. The proposed control 
is implemented for the speed-loop. The simulation results conhrm the effectiveness of the anti-windup controller 
regarding the saturation phenomenon and fast responses. 
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